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ABSTRACT 


Certain ORION missions mav require three axis stabilization. Since ORION’s 
physical size severely limits its onboard fuel storage capability, passive stabilization 
techniques Warrant investigation. This paper shows the development of linearized 
equations of motion and regions of stabilitv with respect to gravitv gradient stabiliza- 
tion. Gravity gradient stabilization by itself provides little yaw restoring torque: there- 
fore, additional torque generating devices are necessary to augment the gravity gradient 
effect. Control moment gvros, reaction wheels, and magnetic torquers will be investi- 
gated as to their suitability for ORION. 
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I. INTRODUCTION 


A. BACKGROUND ON ORION 

The philosophy behind the design and development of the Orion satellite is being 
driven in part by economic considerations. Small, affordable, general purpose satellites 
must be developed without sacrificing reliability. The preliminary groundwork on the 
feasibility of the ORION concept and design is detailed in [Ref. 1}. 

Orion is currently being designed to be launched from a Get Away Special (GAS) 
cannister aboard the Space Shuttle. The GAS cannister physically limits the size of 
ORION to approximately 35” in height and 19” in diameter. Figure 1 shows the Get- 
Away Special Cannister attached to the side of the shuttle bay. 
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Figure 1. N ASA Get-Away-Special (GAS) Cannister [Ref. 1, p. 30] 


The Get-Away-Special (GAS) Cannister 1s designed to launch its satellite 
autonomously on command. A cross sectional view of the (GAS) Cannister is shown in 
Figure 2 on page 2. 

ORION should be able to meet the requirements of standard or typical payloads. 


Some of these requirements are listed by Reference 2 and given in Figure 3 on page 3. 
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Figure 2. Cross section of GAS Cannister [Ref. 1, p. 42] 


1. Background on attitude control 
The success of any satellite is dependent on its attitude control system. To be 
nussion capable, the vehicle must first acquire a desired or specified orientation in space, 
and then maintain this orientaion within mission limits given such external factors as 
magnetic fields, aerodynamic drag, solar pressure, gravitational effects, and other dis- 
turbances. Deviations from desired orientations are detected through the use of sensors. 


Considerations for an attitude control system include: 


Cl) peolnting accuracy 
(2) fuel consumption 
(3) thermal 

(4) power constraints 
(5) satellite lifetime 


(6) response time 
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Figure 3.0 Sumniiary of typical Navy/STP payload requirements 


The weight given to each of these considerations ultimately lies i the specific musston 
requirements designated for each satellite. ‘The direction of this thesis will be to inves- 
tigate the feasibility of achieving sufficient performance utilizing a gravity gradient con- 
trol systent. The standard of +/- 1 degree for each axis will be used. This criteria 1s 
suflicient to enable ORION to satisfy most STP (space test program) mission 
requirements.[Ref. 1, p. 93] | 

‘to achieve stabilization, mamy metliods of generating controlling torques have 
been developed. These methods are divided tnto either active or passive categories. Each 
of these categories has their relative advantages and disadvantages [Figure 4 on page 4 
and ligure 5 on page 4]. 

In this study, the primary focus is on a gravity gradient stabilization augmented 
by other torque generating devices to achieve the desired +/- | degree goal for each axis. 
A discussion of the equations of motion and gravity gradient stabilization theory ts 
presented in Chapter 2. Disturbance torques and other torque generating devices are 
discussed in Chapter 3. Chapter 4 presents the results of stabilization schemes derived 


in Chapters 2 and 3. Observations and conclusions are discussed tn Chapter 5. 
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Figure 4. Passive Stabilization Methods [Ref. 4, p. 19] 
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Fieure 5. Active Stabilization Methods [Ref. 4, p. 19] 


Il. GRAVITY GRADIENT STABILIZATION THEORY 


A. COORDINATE DEFINITION 

Before beginning any discussion on attitude control a means of defining the coordi- 
nate systeni is necessary. The standardized coordinate axes for gravity gradient systems 
are as follows: the x axis points in the direction of satellite motion, the yv axis is normal 
to the orbit plane, and the z axis pomts towards earth's geocenter. These axes corre- 
spond to the roll, pitch and vaw axes, respectively. Figure 6 from Reference 3 defines 
roll (@), pitch (@), and vaw () attitude errors as well as showing the standardized co- 


ordinate system. 
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Figure 6. Coordinate Definition for Gravity Gradient Systems 


B. COORDINATE TRANSFORMATION 


The relationship between two coordinate systems Is essential for attitude control. 


There are inethods available to describe this relationship of which directional cosines, 


quaternions, and Euler angles will be briefly discussed. Figure 7 on page 6 shows various 


methods of coordinate transformations along with their respective advantages and dis- 


advantages. 
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Figure 7, Coordinate Transformation Methods [Ref. 4, p. 412] 
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The directional cosine method involves a 3 x 3 rotation matmx A, where 


where 


Ajj Ay) Ay3 Ax Xey Xe2 


A;; = A>) Ay Ay, —y Oy Yey Yez 


UY 
Ay, Ay) A}; Zex Zey Zw 


A«x=cos@cos 
Xey=cosywsin@singd+sinwycos d 
Xez=sinw sing —cos f sin @ cos @ 
Yex=—sinwcos @ 
Yey=cos@cos py —siny sin@ sing 
Yez=sinwsin@cos¢@+cos w sin 
Zex=sinG 

Zey=—cosésn®g 


ZL *2= COS cose 


(1) 


2. Quaternions 
The quaternion method or as it is sometimes called, the Euler Symmetric Pa- 
rameter method, uses the Euler axis and rotation angle (a single rotation angle « about 
a single axis E) method to describe the relationship between two coordinate systems 


[Figure 8 on page 7]. 


System A 


System B. 





Figure 8. Rotation about the Euler Axis 


The quaternion is a four component vector containing the same information as the Euler 


angle/axis transformation in which one coordinate axis 1s related to another and 1s of the 


form 
Q=L4 Qi, 93k» Gad (3) 
where 
ne 
q, = £, sin > 
eas! 
q2 = Ey sin > 
. oO 
q3 = E, sin > 
44 = COs 2 
2 


The symbols 1, j, and k satisfy the following conditions [Ref. 4, p. 758] 


0 ee 
Jk=—-As=! 
Ai=—ik=j 


and 
2 2) 2 p) 
got+agn+9+q,=1 


When used for control purposes, [Ref 5] defines the first three elements of the 


quaternion (4q,, ¢,. 93) aS the respective roll, pitch, and yaw rotational errors: 


Geran Enitch: eae = 219). 9.930 
= 3 sin, E, sin, E, sin | (7) 
= LOE Gey onaee| 


ow 


The primary advantage of the quaternion method is that the computational time 
can be reduced by more that 40% over the equivalent directional cosine matrix solution 
In Operations that require successive coordinate system rotations [Ref 5]. The transfor- 
mation from coordinate system A to coordinate system B can be described bv the 
quaternion g = [4,/, 9,/. 9k. g, |. and the second transformation from system B to C can 
be described by g’ = [q' |i. g',j. gk. q’,]. The transformation from systenyA to system 


C mav be described by the quaternion'¢ =([¢ Ji, ¢ v0 ong mele 
Q = 4g =(Qit Qa + ak + Gait Vai + 73k + a) (8) 


Muluplving out and substituting in the conditions stated above, the matrix form be- 


cones 
qy G Game 30 4 og leg) 
q'> —-3 Va gq’, Gill @ 
at ae . , (9) 
q 3 G42 —4, Ga 431) % 
Tia} Gd 3 maimed euaae 


Because of the quaternion property of interchangeability, a total of 16 multiplications 
and 12 additions are required to completely define the quaternion describing the trans- 


formation from coordinate svstem A to coordinate system C. To accomplish the same 


transformation using directional cosines would require multiplving two 3x3 matrices, 
involving 27 multiplications and 18 additions. The computational time savings for suc- 
cessive rotations'transformations becomes apparent. The major drawback for the use 
of quaternions is that the numbers of the quaternion do not physically represent the 
transformation from one system to another. 

3. Euler angles method 

The Euler angles utilize three different rotations (angles), defining the orien- 
tation of a body with respect to an inertial reference frame. These angles are defined as 
@ (about the x axis), @ (about the v axis), and & (about the z axis). These axes, by es- 
tablished convention, may also be referred to as | (for x), 2 (for y), and 3 (for z) axes, 
respectively. 

The 321 Euler transformation described below starts with the inertial coordinate 
svstem XYZ. Figure 9 on page 10.a. shows the first rotation about the Z axis through 
an angle y, which produces x,y,Z. Figure 9 on page 10.b. shows the rotation of @ about 
Jy; , producing xj,z,. Figure 9 on page 10.c., showing the rotation of @ about x, produces 
xyz. The rotation described in Figure 9 on page 10.a. about the z axis produces 


Figure 10 on page 10. The transformation from XYZ to x,y,Z 1s 


x, =Xcosyt+ Ysinw 
yp =—-A sinw + ¥cos 


Ld 


aa 
or Written in matrix form 
x) cosw sin OI X 
yy j={-—sinw cosy O}| Y (11) 
vd 0 Oe Liz 





Figure 9. Euler Angles 
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Figure 10. Rotation about the Z axis through the Angle » 


Writing the transformations of Figures 9.b. and 9.c. in matrix form yields 


x cos@ 0 —sin@ ]| x, 
Yi j=} 07 ESO we (12) 
Z, sin@ 0 cos@ jLZ 


and 


x I 0 0 EG 
a Ommcos msi) || 1, (3) 


Z Dieusit D COS.dyai| Z) 


Combining eqs({11-13) into one equation vields the complete 321 coordinate transfor- 


mation from the inertial reference frame to the body reference frame 


x ] 0 0 cos@ 0 —sin@{| cosw sinwW OF} J 
yl|=]0 cos@= sing Ones | 0 —sinw cosy Offi J 
2 O —sing cos¢@]] sn@ 0 cosé 0 0 ILHK 
(14) 
cos @ cos cos @ sin w — sin @ ne 


—cosdsing+singsin@cosw cos@cosw+singdsin@siny sindcosé |i Y 
sngsinw+cos@sin@coswy —singdcosw+cos@siné@sinw cos@cosé }| Z 


This 1s just one of the 12 different Euler coordinate transformations schemes available. 
Figure 11 on page 12 lists all possible Euler angle combinations and their resulting 
transformation matrices. It should be noted that eq(14) does not agree with the 321 
Euler transformation matrix of Figure 11 on page 12. The 321 matrix of Figure 11 on 


page 12 contains a error that interchanges @ with w and w with @. 
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Table of Euler angle combinations [Ref. 4, p. 764] 


C. RIGID BODY EQUATIONS 

The general motion of a particle and its corresponding momentum have been written 
about extensively in various publications. A summary of the relevent derivations and 
equations of motion from [Ref 6] are presented here. 

1. Derivative of a Vector in a Rotating Coordinate System 

The time derivative of a vector in a rotating coordinate system consists of the 

rate of change of the vector relative to the moving axis and the rate of change due to the 
rotation of the axis. Figure 12 shows a vector in a rotating system where XYZ represent 
an inertial system and xyz are the set of axes rotating with angular velocity @ relative 
to XYZ. If unit vectors along the X, Y, and Z axes are i,j, and k, respectively, the vector 


r can be written in the form 


y= xit yj + zk (15) 





Figure 12. Vector in a rotating system [Ref. 6, p. 107] 


The time derivative of r is 
r=xitxit yi+yj +ik+2k (16) 
From [Ref. 7], it has been shown that the time derivatives of the unit vectors are 


i=wxi (17) 


p= wx (18) 
k=wxk (19) 

Defining r,,, as the time rate of change of r relative to xyz, eq(16) becomes 
r="p+wxr (20) 


2. General Particle Motion 


Translational motion of a particle with references to two coordinate systems is 


shown in Figure 13. 





Figure 13. Translational motion of a particle [Ref. 6, p. 107] 


The position vector of a particle relative to coordinate system XYZ 1s 
R=R,+9r (21) 


Where R, is the vector from the axis of the XYZ coordinate frame to the xyz axis frame 
and ris the differential (position) vector in the xyz frame. The velocity of the particle 


is the derivative of eq(21) with respect to time 


V=R=R.+7 (22) 


The velocity at the origin R, = V, combined with eq(20) yields 
V=VitVi,toxr (23) 


3. Momentum of a Rigid Body 
A particle’s momentum is the product of its velocity and mass. The linear mo- 


mentum of the particle of Figure 13 on page 14 1s 


P= | V dm (24) 


where V is the particle’s velocity and m 1s its mass. Recalling that V,,, 1s the time rate 
of change of r relative to the xyz coordinate svstem, and assuming that the xyz axis is 


centered in the rigid body, V,,,=0. Substituting eq(23) into eq(24) with V,,,= 0, 


p= | (Vtaxrdm=mv,40x | ram (25) 


By choosing the center of mass of the rigid body as the center of the rotating coordinate 


system, { r dm=0 by definition. This simplifies eq(25) to 
P=mvlV, (26) 


where P is the linear momentum of the rigid body, m is its mass, and V, is the velocity 
of its center of mass. 
The angular momentum ofa particle dm about point C 1s defined as the moment 


about C of the linear momentum given by 


h,=rxVdm (27) 
Substituting for V yields 
H=| rx Vdm= | rx (V.+0 x) dm (28) 


Rearranging terms 


m 


n 


By choosing the center of mass of the rigid bodv as the origin of the rotating coordinate 


y im) x Vi+ | rxX(@ xan (23) 


m 
system, eq({29) reduces to 


H.=| rx(wxr)dm (30) 


Reducing the angular velocity vector w into its components 
O = Ot + WO + OK (31) 
Combining equations(15), (30), and (31) vields 


o> Uo aoe 
fe yO 1), meee a (32) 


te [1 oO sO tele) Le 


where 


jon -| (s° + 2°) dm (33) 

Ly = | (x° + Za) dm (34) 

I, = | (x? +.y7) dm (35) 
m 


are the moments of inertia about the x,v, and z axes, respectively, and 


Ly = | xy dm (36) 


| | xz dm 
m 


L,= | yzdm 


m 


are the products of inertia. Equation(32) written in matrix form becomes 


fl, x Ly Ly Wx 
Hy | = Ixy ty Aye |] @y 
it, Ly —l,, Lo W, 


4. Equations of motion/Euler’s Moment Equations 


(37) 


(38) 


(39) 


The derivation of the following Euler’s moment equations are summarized from 


[Ref 6]. Recalling that angular momentum /, relative to the point C 1s 

h,=r xX my 
its derivative with respect to time 1s 

h, =rxemv+rrx my 

The moment of the force about a point C bv definition 1s 

Mi=rxF 
Recall Newton's second law describing force: 

F=ma=mvy 


Substituting eqs(22) and (43) into eq(41) yields 


.=rxm(R, +h) +rx F 


Substituting eq(42) into eq(44) and rearranging terms yields the moment equation 


(40) 


(41) 


(42) 


(43) 


(44) 


M,.=H,+ | (R.x F) dm (45) 


m 


where //, is the angular momentum of the rigid body about point C. From eq(20) the 


latter part of eq(45) can be written as 


| R.xrdm=R,x (i (w «| ran) (46) 


m 


Noting that point C is the center of mass f rdm=0 and r,,,=0 for a rigid body, eq(46) 


reduces to 
Mi = Hi. (47) 


The moment of external forces about the center of mass of a rigid body is equal to the 
time rate of change of the angular momentum of the body about the center of mass [Ref. 


7|. In terms of a rotating axis, the moment equation of eq(47) becomes 
M,.= Hoe t+ @ x H. (48) 


Substituting eq(39) in for H, the components of the moment equation can be written as 


M, = H, + w,H, — oH, (49) 
M, = H, + 0H, — oH, (50) 
M, = H, + @,H, — w,H, (51) 


where A/,,, are the moment equations. In this study, an assumption is made that the 
principle axes of the moment of inertia are the rotating axes x,y,z; therefore, the pro- 
ducts of inertia (/,,, /,,, 4,) are equal to zero. The angular momentum components re- 


duce to 


H, = 1,,.0, (52) 


12 EHS Ti 00), (53) 
bE a EO, (54) 


Substituting eqs(52-54) into the moment equations yields Euler's moment equations 


M, = 1p .@y + Wyw(L,, - is) (55) 
ae = Ly @y, + w,o(1,,— 1.) (56) 
M,=1,,0,+ CORO — I.) (57) 


Referring to the coordinate system shown in Figure 6 on page 5, the angular 
velocities in terms of orbital rate w, and attitude error angles pitch(@), roll(@), and yaw 
(w) from [Ref. 6, p. 130] are 


w | 0 — sin 0 db 


x 
@,{=}0 cos@ cos@sin@ 6 
W, 0 —sing cos@cos¢ vy; 
(58) 
cos @ sin w 
—Wy| cosd@cosy+sing sin @ sin 
—singcosy+cos @¢ sin @ sin 
When using small angle approximations, eq(58) reduces to 
Wy db — WoW 
wo, |=| Oe, (59) 
W, W + WoP 
Note that the orbital rate is defined as 
L 
oe (60) 
0 R? 


where yu is the gravitational constant and R is the radius of the earth plus the satellite's 
altitude. For the purpose of station keeping, the small angle approximations are valid. 
Substituting eq(59) into the rigid body moment eqs(55-57) yields the linearized rigid 
body dynamic equations 


Mypg = exh + (hy — Leb — (Ix + Ls — Sy) (61) 


M RB = wv + wo(Iyy ~ Tx) x Og (T x. a I, = Lyd (63) 


D. GRAVITY GRADIENT EFFECT 

Attitude control using gravity gradient stabilization works on the principle that 
given an object with an asymmetrically distributed mass, the object will tend to align it- 
self with its minimum moment of inertia axis along the local vertical towords the earth. 


The force due to gravity on a mass element 1S 


dF =— — R (64) 
where 
R=R)+r (65) 
and 
p= GM, = gR; (66) 


Where G = universal gravity constant, f,= the mass of the earth, g,= gravity acceler- 
ation at the earth’s surface, and R, = the radius of the earth. Reference § shows the 


torque generated about the mass center of the object is 


udm 


Re 





dM,=rxdF=— rx R 


udm 
=— a rx (Ro +7) (67) 





u,am 


mia. 


Cubing eq(65) vields 


3rR, 3r° 3 
?-21+ (68) 
Ro Ro Ro 


20 


Noting that r<R,, the higher order terms drop out. Inversion of eq(68) yields 


3reR 
+=-4+(1-— (69) 
RoR Ri 


Substituting eq(69) into eq(67) and integrating to find the moment yields 


3 x e 
M,=— uel rdmx Ry + | = lal dim (70) 
R 
m 0 m 0 


Given that the center of mass is the origin of the body frame, We dm =(0, eq(70) reduces 


to 





3 
M, = $ | (rx Ro)(r * Ro) dm (71) 
Ro m 


In the reference svstem defined by Figure 14, the K (Z) vector is pointing towards 


the earth; therefore, the vector &, in the reference coordinate frame is 


Ror = RoK (72) 





Figure 14. Gravity gradient torque [Ref. 6, p. 113] 


1] 


Restating the Euler 321 rotation matrix described earlier 


i ] 0 0 cos@ 0 —sin@ |} cosw sinwW O|f J 
Jt=|90 cos@ sing O | 0 —sinw cosw OJ 
k 0 —sing cos¢ {! sin@ O cosé 0 0 LK 
(73) 
cos 6 cos w cos @ sin — sin@ I 


—cos@sing+singsin@cosyw cos@cosw+singsin@siny singcosé@ |i J 
singsinw+cos@sin@cosy —singcosw+cos@sin#ésinw cos¢dcosé }| K 
Transforming R,, into the spacecraft’s fixed coordinates 1, j, and k: 


cos @ cos wy cos 6 sin —sin@ 0 
Ro =| —cos@sing+sing sin@cosw cos@cosw+singsin@siny singcosé@ {| 0 
singsinw+cos@sin@écosyw —singcoswy+cos¢sin#@sinw cos¢dcosé || K 


Meh Teduces (Oo 
Ro, = Ro{ — sin 6i + sin @ cos 67 + cos g cos Bk) (75) 
Substituting eq(75) into eq(71) yields the gravitv gradient moment equations 


(Lo 4,,) Sip Cosi cos’6 
M, = — — (/,, — 1,,) sin 6 cos 6 cos (76) 
2 = (Ty — 1p) Sie eos G-simig, 


“~ 





Replacing 


2a Ue 
Oo = 3 
Ko 





(77) 


where w, 1s the orbital rate, and making the small angle approximations, the gravity 


gradient moment equations become 


Oe = Lo 
M, = 306| (2 — 1,8 (78) 
0 
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Equation(78) clearly shows that by increasing the moments of inertia for /, and /, 
with respect to 7, the gravity gradient moments for roll and pitch will increase. By ex- 
tending a point mass along the z axis, a sizeable increase in the moments of inertia along 
the x and y axes is realized. 


The magnitudes of the gravity gradient torques generated at different altitudes is 
given by Figure 15. 
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Figure 15. Gravity Gradient torque Magnitudes vs Altitude [Ref. 3, p. 58] 





E. DESIGN FOR STABILITY 
1. Stability and Frequency Analysis 
Combining the rigid body dynamic equations (61-63) with the gravity gradient 


equation (78) yields the undamped linear dynamic equations for roll, pitch, and yaw: 
Roll 1,6 + 4wo(l, — 1.) — wo(l, + 1, — LW =0 (79) 


Pitch 1,6 + 3w (I, — 1,0 =0 (80) 


25 


Yaw LW t+ wo(ly— 1) + wll, + 1, - b)b =0 (81) 


Using Reference 9 notations for simplicity, the following variables are defined: 











— J zs 
ye r. (82) 
eae 
p= aT (83) 
Je 
jot 
C= 5 2 (84) 


Substituting eqs($2-84) into eqs(79-S1) yields the following second order differential 


equations 
$ + 4Awod — (1 — A) = 0 (85) 
6+ 3Bw2@ =0 ($6) 
b + Cory +(1- Od =0 (87) 


Taking the Laplace transform of eqs(85-87) allows for frequency analysis. The results 


are put into matrix form: 


Se eV leup 0 —(1—A)wyS |I ¢ 0 
0 S*SaBBa. 0 g|=|0 (88) 
(i= Glaas 0 S*4+Cw, Il y 0 


The determinant of the matrix is the system’s characteristic equation. Using the middle 


row to simplifv calculations, the characteristic equation is found to be 
(S* + (1 +34 + CA)a@gS? + 44 Cag) (S* + 3Bx5) = 0 (89) 


“Solution of the first term on the left vields the natural frequencies of roll and yaw, 
repectivelv; the second term on the left defines the naturai frequency of the pitch axis.” 
[Ref. 9, p. 9] Solving for the first term on the left (roll/yaw stability) bv the quadratic 
formula and dividing through by a} vields 
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2) Saw) ($34 + ACHE A/(1 +34 + ACY ~ 16AC 


Sroll 
ae = (90) 
Cary Wo 5 


For roll/vaw stability the roots of eq(90) must be real and negative; therefore the fol- 


lowing stability conditions exist [Ref. 4 p. 611] 
| 6A AC Se (91) 
AC>0 (92) 


Solving for the second term of Eq(89) and dividing through by @? vields 


Se Ge 
=e Ete) (93) 


e a7 





The condition for stability for the pitch axis 1s 
ASG 2) 


Converting A, B, and C back into their respective moment of inertia relationships de- 
fined by eqs(82-84), combined with the conditions for stability, yields two possible ori- 
entations: | >/,>J/, and /,>J/,>1,. Figure 16 on page 26 shows a summary of the 


stability regions based on moment of inertia ratios. [Ref. 10 p. 204] 
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Figure 16. Gravity Gradient Stability Regions 
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IH. ENVIRONMENTAL DISTURBANCES AND RESTORING TORQUES 


A. ENVIRONMENTAL DISTURBANCES/TORQUES 
The are four major environmental disturbances that affect an orbiting spacecraft: 
aerodynamic, magnetic, gravitv gradient, and solar radiation torque. Figure 17 lists 


these disturbances and their relative regions of dominance. 


SOUNCE OEPENDENCE ON DISTANCE REGION OF SPACE WHERE 


AEROOYNAMIC Pri ALTITUDES BELOW ~ 500 km 


MAGNETIC We ~ 500 km TO ~ 35.000 km: 
GRAVITY GRACIENT We (1 €.. OUT TO ABOUT SYNCHRONOUS ALTITUDE! 


SOLAR RADIATION INDEFENDENT INTERPLANETARY SPACE ABOVE SYNCHRONOUS 
ALTITUOE 


MICNOMETEORITES LANGELY INDEPENDENT; HIGIE CONCENTRATION NORMALLY NEGLIGIBLE; MAY BE IMPORTANT IN 
'N SOME REGIONS OF THE SOLAR SYSTEM SOME SMALL REGIONS TINTERION OF SATURN'S 
RINGS} 


ALTITUDES LISTED AME ONLY REPRESENTATIVE; THE SPECIFIC ALTITUDES AT WHICH VARIOUS TORNOUES DOMINATE ANE HIGHLY 
SPACECRAFT DEPENDENT. 





Figure 17. Environmental Disturbance Torques [Ref. 4, p. 17] 


1. Aerodynamic disturbance 
a. Drag 
“For satellites near the earth the principle non-gravitational force is aero- 
dynamic drag. Aerodynamic drag is a retarding force due to atmospheric friction and is 
in the direction opposite the space vehicles velocity vector.” (Ref. 4, p. 63] 
The drag force is a function of vehicle velocity, air density, and surface area 
of the satellite. The drag force is given by the following equation: 
2 
= pv cya (95) 
2 
where d 1s the drag force, p is the atmospheric density, v is the velocity of the satellite 
along its orbital path, c, is the drag coefficient, and a is the area of the satellite over 
which the drag acts. 
b. Aerodynamic torque 
The collisions of air molecules of the upper atmosphere with the satellite 


surface produces a torque about its center of mass. Aerodynamic disturbances are a 


2) 


function of the satellite’s alutude, velocity, and symmetry. A satellite with a relatively 
low earth orbit may be significantly affected by aerodynamic torques. 

Shown in Figure 18 is a cylindrical satellite with its center of mass (cm) and 
center of pressure (cp). The offset between the center of mass (cm) and the center of 
pressure (cp) is denoted by L,, L,, and L,, respectively. The simplified aerodynamic dis- 


turbance torques from Ref 8 may be computed by the following expression: 
T, = P,1,A, sin & (96) 


Where 7, is the torque due to aerodynamic pressure (ft-lbs), P, is the aerodynamic pres- 
sure (Ibs/f’), L, is the distance between the center of mass and the center of pressure (ft), 
A, 1s the exposed surface area (/?) and a is the angle of attack (radians). For practical 
purposes, the aerodvnamuc torque in the x direction is zero (sina =O when a0). The 


conversion factor from ft-lbs to N-m is approximately 1.356. Figure 19 on page 29 


shows aerodynamic pressure as a function of altitude. 





Figure 18. Aerodynamic effects for cylindrical satellite (Ref. 3, p. 58] 
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NOTE: FROM ATMOSPHERIC DENSITY 
BASED UPON EXTENDED 
ARDC - 1959 MODEL (20) 
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Figure 19. Aerodynamic pressure as a function of altitude [Ref. 8, p. 456] 


The magnitude of the pitch and yaw aerodynamic disturbances for a one 
inch cm-cp offset are shown in Figure 20 on page 30 for various altitudes. The magni- 
tude of these torques decrease significantly with increases in altitude. “Proper control of 
the satellites inertias to enhance gravity gradient torques can make aerodynamic dis- 
turbances essentially insignificant for orbital altitudes above 300 nautical miles (480km). 
On the other hand, in the 100 to 200 nautical mile (160-320km) altitudes aerodynamic 
disturbances are 3 to 4 orders of magnitude greater than the gravity gradient torque from 
a | slug-fr inertia difference and | degree attitude error. In the lower regions, very large 
inertia differences are required for accuracy. [This leads to design considerations] where 
proper configuration can make these aerodynamic torques work as restoring torques 
instead of disturbance torques” [Ref. 3, p. 57]. 
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Figure 20. Aerodynamic torque for a 1 inch offset as a function of altitude. 


2. Rigid Body Solar Pressure Torque 
“The major factors determining the radiation torque on a spacecraft are the in- 
tensity and spectral distribution of the incident radiation, the geometry of the surface 
and its optucal properties, and the orientation of the sun vector relative to the spacecraft” 
BRS emo 7 (0) / 

Solar pressure torques are the disturbances produced by solar radiation pres- 
sure and is a function of the offset of the center of pressure (cp) and the center of mass 
(cm). Referring to Figure 21 on page 31, the center of mass-center of pressure offsets 
are designated as L,, L,, and L,, respectively. The torques generated due to solar pres- 
sure are given by [Ref. 8, p. 454] as: 


T= LA OOo (97) 
T, = 2P,L,A, cos’e (98) 


where 7, = torque due to solar radiation (ft-lbs), P, = radiation pressure (lb/ff?) and is 
~ constant at 9.65 x 10°? lb/fr? for an earth orbiting vehicle, L, = center of mass center 


of pressure offset (ft), A, = surface area of satellite normal to sun(/??), and e = angle 
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of incidence (degrees). Note, equation (97) 1s for an absorbent satellite body and 


equation (98) is for a reflective body. 





Figure 21. Rigid body solar pressure torques [Ref. 3, p. 61] 


B. RESTORING TORQUES 

As seen from chapter 2, the restoring torques of the undamped rigid body combined 
with the gravity gradient effects are considerably smaller for the yaw axis than for either 
the pitch or roll axes. This is due to the necessity of large moments of inertia along the 
pitch and roll axes, as compared to the yaw axis, for gravity gradient stabilization to 
occur. The lack of adequate yaw restoring torque requires that additional yaw restoring 
torque be provided for accurate 3-axis stabilization. A discussion of yaw restoring torque 
methods follows. 

1, Thrusters 

Thrusters are an effective torque generating device; however, for the purpose 

of this thesis, thrusters will not be considered for attitude control. Thrusters are neces- 
sary for orbit maintenance and orbital transfer. Due to the physical size constraints im- 
posed by the GAS cannister, ORION is severely limited in its fuel storage capacity to 
approximately 71 Ibs of hydrazine fuel [Ref 1]. An assumption that the launch platform 
will place ORION in its final orbit is optimistic. Figure 22 on page 32 shows the fuel 
required to go from one orbital altitude to another. Saving the onboard fuel for orbit 
maintenance and transfer and not for attitude control allows more fuel to be devoted 
towards supporting the mission either by allowing greater orbital flexibility and/or a 


longer lifetime. 
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Figure 22. Fuel required for orbit transfer based on a 250I]b satellite. [Ref 2, p. 23] 


2. Momentum Exchange Devices 

The principle behind momentum exchange devices 1s that by changing the an- 
gular momentum of the control device, the angular momentum of the vehicle will be 
changed an equal and opposite amount. Recall that the rate of change of angular mo- 
mentum 1s equal to the torque that ts generated. The two basic types of momentum ex- 
change devices are reaction Wheels and control moment gyros (CMG). The basic 
difference between the two is the method by which the change in angular momentum ts 
accomplished. The reaction wheel axis ts fixed (generally along one of the axes of the 
vehicle) so that the change in angular momentum ts accompltshed by varying the speed 
of rotation of the flywheel. The CMG creates the change tn angular momentum by tilt- 


ing the contant speed flywheel with respect to the vehicle. 
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a. Reaction wheels 
“An axis may be controlled by varying the speed of the reaction wheel in 
response to an attitude error [Ref. 6, p.149]. The momentum of a system with a reaction 


wheel is 
1G SIE Reet (99) 


where H, is the angular momentum of the system and 
Fly x 
jane || Je bre 
H 


is the angluar momentum of the reaction wheel. Aligning the reaction wheels along each 
axis, and substituting eq(99) into eq(48) and adding the gravity gradient moments of 
eq(78) yields 


Inb + 4ere(I, — 1.) — Wolly + I — Lb + (0 - ©) Hys — (W + 08) Hyy + Hux = > Ts 
1,6 + 3000(I,— 1)0 + (— 6 + WW) Hye + (b + 006) Hine + Hy = > Ty (101) 


Lab + wel, — Lb + wo( I, + 1-16 + (6 — wo) + (— 6 + WoW) Myx + Ay: = LE 


“Because H,,,/,,,/1,,, and w, are small, the coupling terms are small. If the coupling 
terms are neglected, the equations of motion about the roll, pitch, and yaw axes become 
independent, and hence, thev can be controlled independently” (Ref: 6, pp. 149-150]. 


Dropping the coupled terms, the motion equations become 


Ib + Anil, — 1) — wo(le + Lp — hb + Hex = ye (102) 
16 + 30%, — 1) + Hy = _T, (103) 
Lb + wi(l, — Lb + wo(lp + —5)b + Hy = )_T; (104) 


oo 


Allowing the momentum of the reaction wheel to be proportional to the attitude error, 


then by differentiation, the reaction wheel control torques are equal to 


Dry = a a) (105) 
Tyrw = — K,(t,8 + 6) (106) 
Taw = ~ K,(t,0 + ¥) (107) 


and the control scheme would be as shown in Figure 23. 


Attitude 
SCNSOTS 


Yaw wheel 





Figure 23. 3 axis reaction wheel control scheme [Ref. 6, p. 149] 


The torque genterated by the reaction wheel is a function of the wheel’s angular velocity 
and its moment of inertia. This means that a small lightweight wheel at a high rate of 
speed generates the same amount of torque that a large wheel at a low rate of speed. 
Since ORION is designed for a relatively short lifespan, wear and tear on spinning parts 
is not an overriding concern. It should be noted that continuous disturbance torques 
that are not cyclical in nature will eventually saturate the reaction wheel, whereby, the 
reaction wheel rotational speed will reach a maximum. Thrusters or some other method 


of momentum removal must be used to reduce wheel speed. Disturbance torques that 
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are cyclical i nature are easily offset bv the reaction wheel as long as the disturbance 
torque is such that it does not saturate the reaction wheel speed. 


3. CMG(s) 


"A control moment gvro’s (either single or double gimbaled) angular mo- 
mentum is due to the rotor which is spinning about the spin axis with a constant 
angular rate. Because the spin axis is gimbaled, a commanded ginbal rotation causes 
the direction of the angular momentum vector to change, thus creating a control 
torque parallel to the output axis. The magnitude of this torque depends on the 
speed of the rotor and the gimbal rotation rate.” [Ref. 4, p. 200] 


A control moment gvro scheme ts shown in Figure 24. 





Figure 24. Twin gyro controller with one drive motor [Ref. 8, p. 415] 


Because of their expense and weight, CMG’s are only used on large spacecraft 


[Ref. 4]. For this reason, control moment gyros will not be addressed. 


4. Niagnetic torquers 

Magnetic torquers may be used to dampen the kinetic energy of the libration 
motion as well as correcting attitude errors in any direction. The magnetic torquers 
consist of three wire coil sets arranged perpendicularly to each other. This configuration 
allows magnetic torque to be generated along any axis or in any direction as a result of 
adding two or three different orthogonal torques. As shown in Figure 25 on page 36, 
when a current 1s applied around the loop of a coil, it produces a magnetic dipole that 
is normal to the plane of the coil and a magnitude which is proportional to the coil’s 
enclosed area and aimpere-turns. 
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n= N/An 


m 





Figure 25. Magnetic Moment Due to a Current Loop [Ref. 4, p. 204] 


The magnetic moment M is defined as 
M=NIAn (108) 


Where N is the the number of turns around the coil, I 1s the current in amperes, A 1s the 
area enclosed by the loop, and 77 is the unit vector normal to the plane of the loop. Be- 
cause of the relationship between N, 1, and A, many different combinations are available 
to produce the desired magnitude of the magnetic moment. Figure 26 on page 37 shows 
three existing spacecraft with various combinations of number of turns and current 


drawn. The torque generated by the magnetic torquer 1s 


Tinag = MX B (109) 


Where M ts the magnetic dipole and B is the earth’s magnetic freld. The earth's magnetic 
field for an equatorial orbit can be approximately modeled as a tilted dipole located at 


78.3° north latitude and 69° west longitude as shown in Figure 27 on page 37. 
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SPACECRAFT REMARKS 


COIL CONSISTS OF 260 TURNS. MAXIMUM 
CURRENT IS 0.6 A AND MAXIMUM POWER 
CONSUMPTION IS 10 W 


COIL CONSISTS OF 360 TURNS. MAXIMUM 
CURRENT 18 0.075 A 


TWO COILS; EACH HAS 500 TURNS. MAX\I- 
MUM POWER CONSUMPTION IS 12 W 





Figure 26. Some existing magnetic moment schemes [Ref. 3, p. 205] 
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Figure 27. Earth’s magnetic field modeled as a tilted dipole [Ref. 6, p. 148] 


The strength of the earth's magnetic field for a circular equatorial orbit in component 


form can be approximated by 


yl 


AM. 


= sin(11.7) sin(281 — lam,) (110) 
Me 
By = ap ca (111) 
A, 
B,= re sin(11.7) cos(281 — lam,) (112) 


Where /am, is measured in degrees of right ascension from Grenich, .M, is the strength 
at the earth’s center which 1s approximately 8 x 10" Wb-m, and R 1s the distance in 
meters from the earth’s center to the orbital altitude. With respect to the satellite, the 
earth's magnetic field 1s changing continuously as the vehicle moves along its orbital 


path. Putting into component form both M and B, eq(109) becomes 


Tinge-x = MB, — M,B, (113) 
Tag-y = MB, — M,B, (114) 
Tmgg-z= MlyB, — M,B, (115) 


It should be noted that the earth’s magnetic field 1s not uniform and varies in 
direction based on the satellite’s position. For accurate attitude control, the vehicle must 
be able to exactly know both its position and the earth’s magnetic field at that position. 
For low earth orbits, the magnetic field is well mapped and may be easily programmed 
into the attitude control computer. For the purpose of this thesis, magnetic field ap- 
proximations will suffice to demonstrate magnetic torquing as a viable attitude control 


capability. 
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IV. CHAPTER 4 


A. GRAVITY GRADIENT STABILIZATION RESULTS 
1. Attitude Control 

There are three general attitude control phases with respect to gravity gradient 
systems: gravity gradient capture, transition period, and steady state. Gravity gradient 
from a tumbling mode occurs when angular velocity slows to one rotation per orbit or 
less, allowing the minimum moment of inertia axis to align with the local vertical. The 
transition period is the time from gravity gradient capture until a steady state has been 
achieved. This period of oscillations will vary depending on the dampening scheme used 
to reduce the amplitude of oscillations. Predominant throughout the satellite’s lifetime, 
the steady state phase is the period in which the satellite’s errors are confined within 
operationally acceptable amplitudes. Maintenance of the desired steady state requires 
minimization of external disturbances and control of those which cannot be eliminated. 

The following group of simulations assumes that a steady state phase has been 
achieved and that the small angle approximations are 1n effect. For each of the simu- 
lations, the following initial conditions are in effect: @, 8, and position error = 5°, and 
d, ’) , and Wy = (), 


2. Various results based on different configurations 
a. Effects of extending booms 
For Figure 28 on page 40, a_10 kg point mass 1s used and the mass of the 
boom is assumed to be negligible. Extension of the boom increases the moment of 1n- 
ertia in the x and vy directions. This increases the restoring torque along the roll and pitch 
axes although it does little to dampen the libration motion. 
b. Gravity gradient effect alone 
The following figures show the gravity gradient effect on a satellite at an 
altitude of 1000 km having moments of inertias as shown. As shown by Figure 29 on 
page 41 and Figure 30 on page 42, stabilization along the roll and pitch axes would take 
a prohibitively long time while the yaw direction would not stabilize. It is clear that 


gravity gradient stabilization alone is incapable of acceptable three axis stabilization. 
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Effect of extending boom on moments of inertia 


Figure 28. 
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Figure 29. Gravity gradient effect only Ix= 15, Iy= 18, Iz=3 
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Figure 30. Gravity gradient effect only Ix=95, ly=98, Iz=3 


c. Effects of Reaction Wheel(s) 
One method of augmenting gravity gradient stabilization is through the use 


of reaction wheels. Teldrix reaction wheel classes and corresponding technical data are 


presented in Figure 31] on page 43. 
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Wheel diameter cm 20 26 35 60°) 


Segesses -ia—— oils a = ee eee. 


Angutar momentuin range Nms 1.8...6.5 14...B80 200... {000 
Max. reaction torque Nin 0.2 0.2 O2 0.3-0.6 
Speed***) min | 6000 6000 6000 
Loss torque at max. speed**) Nin = 0.012 = 0.015 = 0.07 
Power consumption: 


~ steady state (depending Ww se be 2,10 a 10...50 
on speed) 


— max. power rating Ww < <= = 500 


Dimensions: 

— diameter A 350 

— height B 120 | 150 
Weight 5.0...8.0 i ov 2 


Environnrental conditions: 


— operating lemperatiue 

— vibration (sinusoidal) suitable for satellites compatible with laimchers 
— vibration (1andom) such as ARIANE or Space Shuttle 

— linear acceleration 





Figure 31. Teldrix Reaction Wheel Technical Data 


Of the five diameter classes, the 20 cm diameter wheel is the most appropriate of the 
classes listed for ORION. Reaction wheels also come in diameters smaller than 20 cm 
diameter. Bendix Aerospace SP reaction wheel of Figure 32 on page 44 has a diameter 
of © inches. 

As illustrated in Figure 31 and Figure 32 on page 44, power consumption 
values differ with different reaction wheels. For the 20 cm Teldrix reaction wheel, steady 
state requires 2-7 watts with a maximum less than 60 watts. The Bendix SP reaction 
Wheel requires 2.5 watts at steady state with a maximum power consumption of 6.5 
Watts. Due to weight and power consumption values, the Bendix reaction wheel appears 
to be the superior choice of the two. 

The following figures were generated by setting the reaction wheel torques 
equal to eqs(107-109). The initial conditions are @ , 6, and w position errors = 5°, and 
¢=@=W=0. For simulation purposes, the maximum reaction torque value of .2 N-m 
from Figure 31 will be in effect. 

Figures 33-35 have moment of inertias ],=15, J, = 18, and J, =3, with 
K, = K, = 2.0, K, = 1.0, and 1, =1,=1, = 1.0. Figure 33 on page 45 shows ¢ (position 
error), @ (velocity), and the restoring torque (7,,,) generated by the reaction wheel with 
respect to time. Figures 34 and 35 show 6, 6, T.,., and v, wv, Tiwy, Tespectively. With 


the relatively small moments of inertias, stabilization of all axes within the + 1° criteria 
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occurs within 40 seconds. Increasing the moments of incruas to /, = 95, J, = 98 , and 


I, =3 while holding A,, A, at 2.0, A, at 1.0, and 7z,, t, , and t, constant at 1.0, yields 


Figures 36-38. As expected, the @ and @ axes take longer to stabilized (approximately 
160 seconds) than the above case with the smaller moments of inertia. Thé W axis of 
Figure 38 stabilizes much faster than erther the @ or @ axes. This ts due to the relatively 


small moment of inertia (/,) compared to the restoring torque generated by the reaction 
wheel. 


PROGRAM/UNIT NAME : SF Reaction Wheel 


UNIT TYFB NO. (AEF) 3890003 


MOMENTUM (PT-LD-SEC) 0.8 @ 1250 APN 


WRIGHT (LBS. ) 4.9 


SIZ8 DIA. MAX x 3.9" HIGH (excluding aounting 
projection) 


STALL TORQUE (OZ-1N) 2.2 
PEAK FOWER (WATTS) : 6.5 
MAX SPEED (RFR) 1350 


MOTOR EXCITATION : 26V, O-Feak, Squareware, SOO fiz, 2 Phase 


STEADY STATE POWENR( WATTS): 2.5 @ 1350 APM 


TACHOMETER TYPE t P.M. Rotor, 2 Phase Wound Stator, Sinusoidal) 


Output, 2 Cycle/Rev/Phaee, Amplitude varies 
with epeed. 


QUANTITY BUILT ; 82 


COMMENTS AC Squirrel Cage induction Motor 


FLIGHT WESTORY More than 20 of theee units have flown. 


Unite with operating time greater than 5 
yeare are preeentiy Functioning in-orbit. 





Figure 32. Bendix Aerospace Reaction Wheel Technical Data. 
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Figure 33. Gravity gradient effect with 3 RW (X axis) Ix=15, ly= 18, Iz=3 
K=2.0 , T, sel 0 
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Figure 34. Gravity gradient effect with 3 RW (Y axis) Ix=15, ly=18, Iz=3 
K, = 2.0, t, = 1.0 
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Figure 35. Gravity gradient effect with 3 RW (Z axis) Ix=15, Iy=18, Iz=3 
K, = 2.0 , t, = 1.0 
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Figure 36. Gravity gradient effect with 3 RW (X axis) Ix=95, Iy=98, Iz=3 
K,= 2.0, t, = 1.0 
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Figure 37. Gravity gradient effect with 3 RW (¥ axis) Ix=95, ly=98, Iz=3 
K, = 2.0, t, = 1.0 
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Figure 38. Gravity gradient effect with 3 RW (Z axis) Ix=95, ly=98, Iz=3 
K,=1.0,7,= 1.0 
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Figures 39-41 have moments of inertia of J, = 15, 4, = 18, and 7, =3 with 
K, = K, =.2, K,=.1, and t, =1, =t,= 1.0. The reduction of A,, A, and K, reduces the 
reaction wheel torque generated. 

Figure 39 shows (position), é(velocity), and reaction wheel uenmefue { 
with respect to time. Figures 40 and 41 show @, é, tPeeancneay, | TeESDeClively, 
Stabilization, within the + 1° criteria, occurs for all axes in under 200 seconds. Com- 
parison of the magnitudes of the torque values of Figures 33-35 and Figures 39-41 for 
the same moment of inertia values shows a decrease of roughly a factor of five. 

Increasing the moments of inertia to /, = 95, /, = 98 and J, = 3, while hold- 
ing A, = K,=.2, K,=.1, and t,=7t,=t, = 1.0, produces Figures 42-44. This last group 
of reaction wheel figures shows that +650 seconds are required to dampen the system 
to the + 1° stabilization criteria. The tradeoffs between the different reaction wheel 
combinations are essentially time response against torque/power. The faster the desired 
response, the greater the magnitude of the restoring torque, and ultimately, the more 


power that will be consumed in achieving this response. 
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Figure 39. Gravity gradient effect with 3 RW (X axis) Ix=15, ly=18, Iz=3 
Ke = 2. = 10 
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Figure 40. Gravity gradient effect with 3 RW (Y axis) Ix=15, Iy=18, Iz=3 
K,=.2,t1,= 1.0 
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Figure 41. Gravity gradient effect with 3 RW (Z axis) Ix=15, ly=18, Iz=3 
K,=.1,7,=1.0 
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Figure 42. Gravity gradient effect with 3 RW (X axis) Ix=95, Iy=98, Iz=3 
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Figure 43. Gravity gradient effect with 3 (Y axis) RW Ix=95, Iy=98, Iz=3 
K,=.2,7,= 1.0 
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Figure 44. Gravity gradient effect with 3 RW (Z axis) Ix=95, ly=98, Iz=3 
ko ——lapaoe— 1.0 
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d. Effects of Magnetic torquers 

For the magnetic torquing simulations, the initial conditions are ¢, 0, and 
W position errors = 5°, and ¢. 6. and w = 0. The initial restoring torque is provided 
by the gravity gradient effect combined with the rigid body dynamics. For the following 
simulations, one of the two conditions must exist prior to energizing the magnetic coils 
to generate magnetic torque: 

1) position greater than .5° and velocity greater than 0, or 

2) position less than —.5° and velocity less than 0. For the following figures, 
the magnetic moment varies from 9.0427 to .5 Wb-m. At the maximum value of 
9.0427W-m, N= 400 turns, [=.5 amp, and A= zr? where ris the radius at .12m. Subse- 
quent values are obtained by reducing the current drawn. 

The magnetic field in component form for an equatorial orbit can be ap- 
proximated by eqs(110-112). Since the earth’s magnetic field is highly dependent on the 
satellite's position. ORION must have the capability to determine its position in order 
to effectively utilize magnetic torquing. For simulation purposes, the satellite 1s assumed 
to be directly over Greenwich at time t=0. The magnetic coils are turned on and off at 
the specified strength as needed. Varving the strengths of the magnetic moments gener- 
ated would enhance performance, but is bevond the scope of this thesis. The assumption 
that magnetic torque may be generated in any direction due to vector addition is not 
alwavs valid: magnetic torque may not be generated in the direction of the magnetic 
field. 

Figures 45-47 utilized magnetic moments of m, =m, =m, = 9.0427 with 
moments of inertias /, = 15, J, = 18, and /,=3. For Figure 45 on page 60, the gravity 
gradient restoring torque takes approximately 800 seconds before meeting a condition 
necessary for energizing the magnetic coils. Due to the small moments of inertia, the @ 
restoring torque is needed for only a short time. Figures 46 and 47 show the @ and w 
axes essentially stabihzed within the + 1° criteria after the the initial magnetic torque 
pulse. For the @ axis, this pulse occurs at the 1200 seconds and lasts for approximately 
300 seconds. The @ axis has a lower frequency response and, therefore, takes longer to 
to meet the condition necessary to energize the magnetic coils. Due to its relatively small 
moment of inertia, the ¥ axis chatters slightly before stabilizing. A smaller magnitude 
of restoring torque eliminates this problem. 

Increasing the moments of inertia to J, = 95, J, = 98, and J, =3 while leav- 
ing the magnetic moments m,, m,, and m, constant at 9.0427wb-m, produces Figures 


48-50. As expected, more restoring torque is needed to stabilized the @ and @ axes forc- 
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ing the magnetic coils to remain energized for a longer period of time. The time needed 
for all axes to be within the + 1° criteria 1s limited by the @ axis. Stabilization occurs 
within 2800 seconds while for the case with the smaller moments of inertias, stabilization 


occurs at approximately 2100 seconds. 
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Figure 45. Gravity gradient with magnetic torquing (X axis): Magnetic 
moment = 9.0427W-m Ix= 15, ly= 18, Iz= 3. 
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Figure 46. Gravity gradient with magnetic torquing (Y axis): Magnetic 
moment = 9.0427W-m [x= 15, ly= 18, Iz= 3. 
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Figure 47, Gravity gradient with magnetic torquing (Z axis): Magnetic 
moment = 9.0427W-m [x= 15, ly= 18, [z= 3. 
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gradient with magnetic torquing (X axis): Magnetic 
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Figure 49, Gravity gradient with magnetic torquing (Y axis): Magnetic 


moment = 9.0427W-m Ix= 95, ly= 98, Iz= 3. 
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Figure 50. Gravity gradient with magnetic torquing (Z axis): Magnetic 
moment = 9.0427W-m Ix=95, ly = 98, Iz= 3. 
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Reducing the strength of the magnetic moments reduces the torque that is 
generated and, thereby, slows the response time to stabilization. For Figures 5l-o0ntne 
magnetic moments m, = m1, =m, = 1.0. Figures 51-53 have moment of inertias J, = 15, 
I,=18 , and /,= 3. Stabilization for the @ axis occurs within approximately 3600 sec- 
onds. Increasing the moments of ineruas to J, = 95, /, = 98, and /, = 3 produces Figures 
54-56. Stabilization within the + 1° criteria for the @ axis occurs at approximately 11000 
seconds. The length of time that the magnetic torgers are energized increases as the 
strength of the moment decreases. Figures 57-62 show the effects of magnetic torquing 
using a maximum magnetic moment of .SWb-m for each coil. Figures 57-59 show that 
for low moments of inertia (/, = 15, J, = 18, and /, = 3), stabilization occurs within the 
+ 1° within 4800 seconds. Increasing the moments of inertia to 1, =95, J, = 98, and 
I,=3 holding the maximum magnetic moment =.5Wb-m, produces Figures 60-62. In 


Figure 61, the @ axis may take a prohibitively long ume. 
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Figure 51. Gravity gradient with magnetic torquing (X axis): Magnetic 
moment = 1.OW-m Ix= 15, ly= 18, Iz=3. 
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Figure 53. Gravity gradient with magnetic torquing (Z axis): Magnetic 
moment = 1.0W-m [x= 15, ly= 18, [z= 3. 
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Figure 55. Gravity gradient with magnetic torquing (Y axis): Magnetic 
moment = 1.0W-m Ix=95, ly= 98, Iz=3. 
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Figure 56. Gravity gradient with magnetic torquing (Z axis): Magnetic 


moment= 1.OW-m Ix=95, Iy= 98, Iz= 3. 
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Figure 538. Gravity gradient with magnetic torquing (Y axis): 
moment = .S5O0W-m [x= 15, ly= 18, Iz=3. 
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Figure 59. Gravity gradient with magnetic torquing (Z axis): Magnetic 


moment= .50W-m Ix= 15, ly= 18, Iz=3. 
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Figure 60. Gravity gradient with magnetic torquing (X AXIS): Magnetic 


moment = .50W-m 1x=95, ly= 98, Iz= 3. 
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Figure 62. Gravity gradient with magnetic torquing (Z AXIS): Magnetic 
moment= .50W-m Ix=95, ly = 98, Iz= 3. 
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3. Gravity Gradient capture (from tumbling mode) 

Increasing the inertia will reduce the tumbling rate by conforming to the con- 
servation of angular momentum, /,@, = /,w,. As the inertia increases, a corresponding 
decrease in angular rate must be realized. This 1s useful in obtaining the initial gravity 
gradient capture from a tumbling mode. Figure 63 on page 80 shows that for a slow 
tumbling satellite, the change in the moment of inertia by extending the boom may be 
sufficient to allow gravity gradient capture. The satellite is initially tumbling at a rate 
of 5 rotations per orbit about the pitch axis. The extension of the boom increases the 
moment of inertia about the y axis by a factor of 6; therefore, by the conservation of 
angular momentum, the angular rate must decrease by a factor of 6. This reduces the 
rotations per orbit to approximately 5/6 per orbit well below the one rotation per orbit 
rate necessary to acheive gravity gradient capture. Figure 63a shows the initial tumbling 
rate of 5 rotations per orbit. Figure 63b shows the corresponding decrease in the tum- 
bling rate as the moment of inertia 1s increased. Figure 63c shows the increase in mo- 
ment of inertia as the boom 1s extended with respect to time. 

When an initial tumbling rate is higher than the rate that can be captured solelv 
through the extension of the boom, magnetic torquing may aid in gravity gradient cap- 
ture. By creating torque in the direction opposite to that of the angular velocity, mag- 


netic torquing can be used to slow the rate of the tumbling vehicle. 
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Figure 63. Gravity gradient capture by extension of boom 


80 


Vo CHAPRRR 5 


A. CONCLUSIONS 
1. Gravity gradient stabilization (only) 

Gravity gradient stabilization alone is not a viable attitude control method for 
ORION. This is due to gravity gradient effect providing little yaw restoring torque. In 
addition, gravity gradient stabilization provides little dampening effect on the roll and 
pitch axes. When augmented by other control schemes, gravitv gradient stabilization 
provides sufficient passive restoring torque to warrant its use. 

2. Gravity gradient with 3 reaction wheels 
a. Advantages 
Gravity gradient stabilization augmented by three reaction wheels, each 
aligned along the body’s axes, provides excellent attitude control results both in time 
response and accuracy. 
b. Disadvantages 
Power consumption may be prohibitively high for this method. At steady 
State, each reaction may use up to 6 watts while drawing considerably more during the 
acceleration or deceleration phase. If a steady state drain of 18 watts (6 for each re- 
action wheel) 1s acceptable, then this method should be considered. The projected end- 
of-life power budget for ORION is 60 watts [Ref. 1, p. 132]. Using this figure as an 
estimate, a steady state draw of 18 watts leaves 42 watts available for mission support. 
From Figure 3 on page 3, an average of 34 watts is required for mission support. The 
size and weight of three reaction wheels may be prohibitively high as well. Another 
disadvantage is that thrusters must be used to remove momentum when the reaction 
wheel(s) become saturated. 
3. Control moment gyro (CMG) 
a. Advantages 
Although not discussed in detail, control moment gyros are capable of 
generating large restoring torques. 
b. Disadvantages 
The weight requirements for a control moment gyro are prohibitively high 


when compared to reaction wheels when relatively small torques are invoved. As a small 
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multi-purpose satellite, ORION does not need the magnitude of restoring torque neces- 
sary to warrant the use of control moment gvros. 
4. Gravity gradient with three magnetic torquers 
a. Advantages 
Three othogonally oriented magnetic torquers have the advantage of gen- 
erating restoring torque in almost any direction based on vector addition. In this man- 
ner, the corrective torque may be produced to counter any disturbance from any 
direction as long as it is small enough not to tumble the vehicle. Another advantage for 
magnetic torquing 1s that the torquers only draw power while providing restoring torque. 
This amounts to a considerable power savings When compared to the reaction wheel. 
The power drawn by the magnetic torquers during operations 1s based on the amperes 
used. This allows for flexibilitv of design since doubling the number of turns and halving 
the current drawn still provides the same magnitude of the magnetic moment. 
b. Disadvantages 
Magnetic torquers are altitude dependent. The magnitude of the earth's 
magnetic field is a function of ry At high altitudes, magnetic torquing Is not effective. 
In addition, magnetic torquing is not always available in all directions. Recalling that 
magnetic torque is 7,,, = .\/ x B, magnetic torque may not be generated in the direction 
of the B field due to vector cross product relationships. Due to the continuously chang- 
ing B field with respect to the satellite, this is a minor disadvantage. The desired restoring 
torque can soon be generated as soon as the B field direction changes. 
5. Thrusters 
a. Advantages 
At high altitudes, thrusters and momentum exchange devices are the only 
control methods available for accurate stabilization. In addition, thrusters are the most 
effective method for stabilizing a vehicle with a high tumbling rate. 
b. Disadvantages 
Thrusters are necessarv for satellites; however, anv fuel used for attitude 
control or momentum removal is less fuel available for orbit transfer and maintainance. 
For a small, lightweight satellite with limited fuel storage capacity, fuel resources must 
be economized. 
6. Overview 
In the final analysis, the attitude control system employed will be based on 


specific mission requirements. Only when the ORION mission is Known can proper 
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weights be given to each of the attitude control system design factors, most noticeably, 


power consumption, pointing accuracy, altitude, and desired lifetime. 
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APPENDIX GRADIENT STABILIZATION PROGRAM 


TITLE 3 AXIS STABILIZATION WITH GRAVITY GRADIENT EFFECT 


PTT TET TET TETETETS TEE TS ETE TEES BEET TE RETEST ETE TSE EEREEIERREERIER EDEL ESE 


7 DEFINITION Oise TERMS 


* ALT=ALTITUDE OF SATELLITE (METERS) 
* R1=RADIUS OF EARTH (METERS) 
** R=RADIUS OF EARTH + SATELLITE ALT (METERS) 
* RO=DENSITY OF AIR AT THE ORBITAL ALTITUDE 
* WO=ANGULAR RATE OF SATELLITE WITH RESPECT TO EARTH 
*  D=DRAG FORCE=1. 7E+14*(RO/R) 
* PSI,THETA,PHI ARE ANGLES ABOUT THE X,Y,Z, AXES (RADIANS) 
* PSIC,THETAC,PHIC ARE ANGLES ABOUT THE X,Y,Z AXES (CDEGEREBo 
* PSID,TD,PHID ARE ANGULAR VELOCITIES ABOUT THE X,Y,Z AXES (RADS7Shae 
* (B1),BC2),B(3}) ARE THE ANGULAR ACCELERATION ABOUT THE 
* X,Y,Z2 AXES (CRADS/SEC**2 ) 
* IX,IY,I1Z ARE THE MOMENTS OF INERTIA ABOUT THE X,Y,Z AXES (KG-M**2) 
* TX,TY,TZ ARE THE SUMMATION OF THE DISTURBANCE TORQUES AND 
¥ RESTORING TORQUES EXCLUDING GRAVITY GRADIENT EFFECTS 
* TMX, IMY,INZ ARE MAGNETIC TORQUES ABOUT THE. x12 752 53 
* TRWK,TRWY,TRWZ ARE REACTION WHEEL TORQUES ABOUT THE X,Y,Z AXES 
*  TSX,TSY,TSZ ARE SOLAR PRESSURE TORQUES ABOUT THE Xoy 34a 
* TAX,TAY,TAZ ARE AERODYNAMIC TORQUES ABOUT THE X,Y,Z AXES 
* LX,LY,LZ IS THE DISTANCE BETWEEN THE CENTER OF (PRESSURE ane 
7 CENTER JOP MAS Sects 
* MX ,MY,MZ=SPACECRAFT'S GENERATED MAGNETIC MOMENT WITH RESPECT TO 
¥ THE X,Y,Z AXES (WB-M} M=NIA WHERE M=MAGNETIC MOMENT, 
We N=} OF TURNS, I=CURRENT (AMPS), AND A=AREA (M***2). FOR 
¥ THE MAXIMUM VALUE OF M=9.0427, N=400,I=. 5A,RADIUS=. 12M. 
“= MG=STRENGTH OF EARTH'S MAGNETIC DIPOLE=8.0E15 WB-M AT R=0 
* BX,BY,BZ=EARTH' S MAG FIELD WITH RESPECT TO X,Y,Z(WB/M***2) 
¥ LAMS=RIGHT ends dee aioe Gidea ines @ TIME=0, IN DEGREES( LONGITUDE) 
Yerteleds ves stestes feveste ete slevedesleste s'e s'e s'es nto evesles'e eslestc sles'e slorlestos owe evtesedevesleve Jenleslenietee ge Tatras s'e sesleveveve 
ARRAY B(3) 
CONST ALT=1000E+03, RO=7. OE-10,R1=6378. E+03,... 
[X=95., 1Y=968., 12=03.00, TANso-0-0. 
MX=9.0478, MY=9.0427, MZ=9.0427, MG=8. OE+15,... 
KX=2., KY=2., KZ=1., TAUI=1., TAUZ=10 se See 
TSX=0. , TSY=0., TSZ=0. , TAX=0, , TAY=O0- ee —O- 
* 
dededetesesetevekdesetetetesetevetetetetetetetetedehedetetetetedede tevete detec Rae he teak 1c sete eto ie is ae ae Teak oe nee 
* PHI , THETA, AND PSI ARE INIT 5 DEGREES = .08725 RADIANS OF ERROR* 


io oe ef o omle elewdeutes eewewe een ee owowsesuwou'e ale wen’ ou'e = oe Pd - io o ate oe eowou' es owe 
vs Seveveds seve Jeveys vev's e ve qe =e rts we =e vs ws ve re “¢ we ws vs =" ss ws ms wie ¥ =‘ se se Ae mses nevtew's desles's ves vevess ve qe Yevs seve oe ss vesevevevs ws veveveve 


INCON ITHETA=. 08725 , IPHI=. 08725 ,IPSI=. 08725 , ILAMS=0 , PHIC=5. , 


THETAC=5. ,PSIC=5. 
fe 


DERIVATIVE 


=R1+ALT 
D=1. 7E+14*(RO/R) 
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WO=SQRT( 3. 987E+14/(R**3)) 


NOSORT 
LAMS=LAMS0+WO*TIME*(57. 3) 

10 itt CLAMS. GE. 35600) THIEN 
LAMS=LAMS -360. 


co TO) 10 

ENDIF 
* 
* 
dedededededeacdedkedcskakvdcacdedededededescdedk tedevede deve dedededede deve desde sevens tok seve dosed tek deve dex 
* THIS SECTION ADDS THE REACTION WHEEL EQUATIONS 
Oe A IA Se ee ie ae oe be ae oie ie ie ee eo oe ie Be EL ERR DIR Re ta eR ie eed 
* TRWX=-KX*( TAU1*PHID + PHI) 
* TRWY=-KY*( TAU2*TD + THETA) 
* TRWZ=-KZ*( TAU3*PSID + PSI) 


Tedkdeakedk deve dsaedewt F NT) REACTION WHEE L228 28% ae 2 ve ae ae dere eve Teak ae He sek Heke ak Ie Fe 
sedeseaedcseacseskolcacdedeseskvokvdcacdevcacdeseckakdessdeacsssesedcacakcseseslscdestacscseskoteskdkdkeveseskaksedesk 


* THIS SECTION ADDS THE MAGNETIC TORQUE EQUATIONS * 
* AND SOLVES AN EXAMPLE FOR THE EARTH'S MAGNETIC * 
* FIELD (EQUATORIAL ORBIT ASSUMED) * 


2 ik oer he oe ae a bir bt ik bee bi oe bie oir ie ir oir ie ok bie bk oie ok bie) teat ee ie oe ok $e Oe ie a ke oe ok ke oie te ne ae 
BX=(MG/(R**3))*SINC. 2042)*SIN(C (CC LAMS*PI)/180)-5. 0789) 
BY=-( MG/(R**3))*COSC(. 2042) 
BZ=-(MG/CR**3) )*SINC. 2042)*COS( ( (LAMS*PI)/180)-5. 0789) 


rly 
TMAGX=( MZ*BY) -(MY*BZ) 
TMAGY=( MX**BZ) -(MZ**BX) 
TMAGZ=( MY*BX) -(MX*BY) 


Sevesetescscdevevesesesevedesesesedesesesdevete dese de sede seek 


* CONDITIONS FOR TORQUE GENERATION MET? 


Seicdesesievesedesevessdevedesesevedsdetsvsvededkdelsvededes's 


IF(PHIC.GT..5) THEN 
IF(PHID. GT. 0) THEN 
TMX=-( ABS( TMAGX) ) 
ELSE 
TMX=0 
ENDIF 

ENDIF 

IF(PHIC. LT. -.5) THEN 
IF(PHID. LT. 0) THEN 
TMX=ABS( TMAGX) 

ELSE 
TMX=0 
ENDIF 

ENDIF 

IF(THETAC. GT..5) THEN 
IF(TD.GT.0) THEN 
TMY=~(ABS( TMAGY) ) 
ELSE 
TMY=0 
ENDIF 

ENDIF 

IF(THETAC. LT. -.5) THEN 
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IF(TD. LT. 0) THEN 
TMY=ABS( TMAGY) 
ELSE 
TMY=0 
ENDIF 

ENDIF 

[F(PSICHGIey one HEN 
IF(PSID.GT.0) THEN 
TMZ=-(ABS(TMAGZ) ) 
ELSE 
TMZ=0 
ENDIF 

ENDIF 

TF CPSlO spines nen 
IF(PSID. LT.0) THEN 
TMZ=ABS(TMAGZ) 
ELSE 
TMZ=0 
ENDIF 


ENDIF 
Hawes severedesese ves sewer wee KENT MAGNETIC TORQING SECTION tot dedvede ae ied desedevededed dese 


TX=TMX+TRWA+TSX+TAX 
TY=TMY+TRWY+TSY+TAY 
TZ=TMZ+TRWZ+TSZ+TAZ 


Weaker OR GRAVITY GRADIENT EFFECT ONLY SET TX=TY=T Z=O2' ttc cevede deat vevtede secede 


BC 1)=-4*(WO**2)*( CTY-IX)/1Z)*PHI+WO*( CIX+1IZ-1IY) /IX)*PSID+(TX/I1X) 
B( 2)=-3*%(WO**2)*( (IX-1Z)/IY)*THETA+TY/IY 
B(3)=-(WO**2)*( (IY-1X) /1Z)*PSI-wO*(( IX+1Z-IY)/1Z)*PHID+(TZ/1Z) 
ry 
PHID=INTGRL(0. ,B(1)) 
TD=INTGRL(CO. ,B(2)) 
PSID=INTGRL(0. ,B(3)) 
* 


ve 
CPHI=INTGRL(CO. , PHID) 
CTHETA=INTGRL(O. , TD) 
CPSI=INTGRL(O. , PSID) 


PHI={1PHI+CPHI 

THETA=ITHETA+CTHETA 

Pob=lPSi+Crsi 
detevededevedevevevedesesesodesedeseds sede tevedede lode de dede code fete sede ese deveT RK FIT OTC TTC Te eee 
*PHI,THETA,PSI ARE = TO INITIAL CONDITIONS AS ANNOTATED WITH AN "I" 
*PRECEDING VARIABLE+CALCULATED VALUES AS ANNOTATED WITH THE "C" IN FRONT 
dadvdevedededcdsdeds dededededede ve de deve vedo dk dese dedevededevededede deve vevedevedeveve de deve Fede dete seve Fete te dete se sede seve sete vee sede seus a 

PHIC=PHI*57. 3 

THETAC=THETA*57. 3 

PSIC=PSI*57. 3 


ri 
Tek deddedededesevededededededdevedevevesevedededtededevdedvvdevestededeseseslevdescseseotesesededeseslestesle se seve sleve skate se se dese ake ve se ae 
* USE FOLLOWING IF GRAVITY GRADIENT EFFECT ONLY ¥ 
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Jeseserevicdededevedededesevededesesedeskedeskesesedevededesdcdeddevedtdcdededdededededs sede vetevedetevedvesvededesededvevededeveveveveve 


*ONTRL FINTIM=20000 ,DELT=200. 

“RINT 200., PHIC, THETAC, PSIC, PHIDC, TDC, PSIDC 

*AVE B00]. PHIC THETAC] PSlGeee ire, TDC, PSIDC 

*RAPH (G1,DE=TEK618) TIME( UN=SEC) , PHIC( UN=DEGREES ) 

*ABEL (G1,DE=TEK618) PHICDEG) VS TIME (GG@1000K) IX=95 IY=98 IZ=3.00 
*RAPH (G2,DE=TEK618) TIMECUN=SEC) , THETAC( UN=DEGREES ) 

*ABEL (G2,DE=TEK618) THETACDEG) VS TIME (GG@1000K) IX=95 IY=98 IZ=3. 00 
*RAPH (G3,DE=TEK618) TIME(UN=SEC) , PSIC( UN=DEGREES ) 

*ABEL (G3,DE=TEK618) PSICDEG) VS TIME (GG@1000K) IX=95 IY=98 IZ=3.00 
*RAPH (G4,DE=TEK618) TIME( UN=SEC) , PHIDC( UN='DEG/SEC' ) 

*ABEL (G4,DE=TEK618) PHID VS TIME (GG@1000K) IX=95 IY=98 IZ=3. 00 
*RAPH (G5 ,DE=TEK618) TIME(UN=SEC) , TDC(UN='DEG/SEC' ) 

*ABEL (G5,DE=TEK618) TDC VS TIME (GG@1000K) IX=95 IY=98 IZ=3. 00 
*RAPH (G6,DE=TEK618) TIME(UN=SEC) ,PSIDC( UN='DEG/SEC' ) 

*ABEL (G6, DE=TEK6 18 ) PSIDC VS TIME (GG@1000K) IX=95 IY=98 IZ=3. 00 
veveredededededeveseseseKEND GRAVITY GRADIENT ONLY GRAPHS 3 Xaevtdsdedeveuededsdedevedevedeue deve 

% 

Dae ae as Wea es I FC 20 oe BC OC TT OC IC DE ee IC I IE Ie TENE TC IE OC IE Ie Oe ae IE OK eae Ie eae IC Ie IE Oe Ie IE Ie IC Ie Ie IS Oe TE Oe Mee NS 

* USE THE FOLLOWING IF REACTION WHEELS ARE USED * 

Bee Fee Ie I 90 96 1S FS IC ICT NE IS TE TENS FETC IE TE TOTTI IE TE TE IE TE We ETE IE IE IE TENE I IE IE ISDE IME DENE DETE DE DEE BENE INE 

*ONTRL FINTIM=400,DELTS=5. 

*RINT 100., PHIC, THETAC, PSIC, TX, TY, TZ 

“AVE 5., PHIC, THETAC, PSIC, TX, TY, TZ, PHID, TD, PSID 

*RAPH (G1,DE=TEK618) TIME(UN=SEC) , PHIC( UN=DEGREES ) 

*ABEL (G1,DE=TEK618) PHI(DEG) VS TIME (GG@1000K) IX=95 IY=98 IZ=3.00 
*RAPH (G2,DE=TEK618) TIME(UN=SEC) , THETAC( UN=DEGREES) 

*ABEL (G2,DE=TEK618) THETAC( DEG) VS TIME (GG@1000K) IX=95 IY=98 IZ=3.00 
*RAPH (G3,DE=TEK618) TIME( UN=SEC) , PSIC( UN=DEGREES ) 

*ABEL (G3,DE=TEK618) PSI(DEG) VS TIME (GG@1000K) IXK=95 [Y=98 IZ=3. 00 
*RAPH (G4,DE=TEK618) TIMECUN=SEC ) , TX( UN=NM ) 

*ABEL (G4,DE=TEK618) TXRW VS TIME (GG@1000K) IX=95 IY=98 IZ=3.00 
*RAPH (G5,DE=TEK618) TIME(UN=SEC) , TY(C UN=NM) 

*ABEL (G5,DE=TEK618) TYRW VS TIME (GG@1000K) IX=95 IY=98 IZ=3. 00 
*RAPH (G6,DE=TEK618) TIMECUN=SEC) , TZ( UN=NM) 

*ABEL (G6,DE=TEK618) TZRW VS TIME (GG@1000K) IX=95 IY=98 IZ=3. 00 
**RAPH (G7,DE=TEK618) TIME(UN=SEC) , PHID(UN='RAD/SEC' ) 

*ABEL (G7,DE=TEK618) PHID VS TIME (GG@1000K) IX=95 IY=98 IZ=3. 00 
*RAPH (G8,DE=TEK618) TIME(UN=SEC),TD(UN='RAD/SEC' ) 

“ABEL (G8,DE=TEK618) TD VS TIME (GG@1000K) IX=95 IY=98 IZ=3.00 

*RAPH (G9,DE=TEK618) TIME(UN=SEC) , PSID(UN='RAD/SEC' ) 

“ABEL (G9, eee) PSID VS TIME (GG@1000K) IX=95 IY=98 IZ=3.00 
Wresodedededs Hee seTe eves AK EE NT REACTION WHEEL GRAPHS vt 2veadevededeseveds Teaedesesesese dose 
we ve 
Jered I sesesvetesetede doves dots seteve sede dove deve sede se Ke eT Tete te KR ete deve te Fete se Fe dedede Fe eek 
* USE FOLLOWING GRAPHS IF MAGNETIC TORQUING IS UTILIZED * 
Wea everesesesese sere sel seve se Tove ee Fe se Fe Fe Fe MeN FE He Fee Te Tes We Fe Ve He FE Te Te Te Fe Fe Fe Ne Vee Fe Te Fe Fe Te TENT ARN Tevete eI 
*ONTRL FINTIM=4000. ,DELT=5. 

*RINT 100., PHIC, THETAC, PSIC, WO,B(1),B(2),B(3),LAMS 

*AVE 5., PHIC, THETAC, PSIC, LAMS, TX,TY,TZ,BX,BY,BZ,PHID,TD,PSID 
**RAPH (G1,DE=TEK618) TIME( UN=SEC) , PHIC( UN=DEGREES ) 

*ABEL (G1,DE=TEK618) PHI(DEG) VS TIME (GG@1000K) IX=95 IY=98 IZ=3. 00 
*RAPH (G2,DE=TEK618) TIME(UN=SEC) , THETAC( UN=DEGREES) 

*ABEL (G2,DE=TEK618) THETACDEG) VS TIME (GG@1000K) IX=95 IY=98 IZ=3. 00 
*RAPH (G3,DE=TEK618) TIME(UN=SEC) , PSIC( UN=DEGREES ) 

*ABEL (G3,DE=TEK618) PSI(DEG) VS TIME (GG@1000K) IX=95 IY=98 IZ=3. 00 
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*RAPH 
*ABEL 
*RAPH 
*ABEL 
*RAPH 
*ABEL 
*RAPH 
*ABEL 
*RAPH 
*ABEL 
*RAPH 
*ABEL 
*RAPH 
*ABEL 
*RAPH 
* ABEL 
*RAPH 
*ABEL 
**RAPH 
* ABEL 
END 

STOP 


(G4, DE=TEK618) 
(G4 ,DE=TEK618) 
(G5 ,DE=TEK618) 
(G5 , DE=TEK618) 
(G6 ,DE=TEK618) 
(G6 , DE=TEK618) 
(G7 , DE=TEK618) 
(G7, fh=—TEKele) 
(G8 , DE=TEK618) 
(G8 , DE=TEK618) 
(G9, DE=TEK618) 
(G9 ,DE=TEK618) 


(G10 ,DE=TEK618) 
(G10 ,DE=TEK618) 
(G11,DE=TEK618) 
(G11,DE=TEK618) 
(G12, DE=TEK618) 
(G12 ,DE=TEK618) 
(G13 ,DE=TEK618) 
(G13 ,DE=TEK618) 


TIME(UN=SEC) , LAMS( UN=DEGREES ) 
LAMS(DEG) VS TIME (GG@1000K) IX=95 IY=98 I2=3.00 
TIME(UN=SEC) , TX(UN=N-M) 
TX VS TIME (GG@1000K) 
TIME( UN=SEC) , TY(UN=N-M) 
TY VS TIME (GG@1000K) 
TIME(UN=SEC) , TZ(UN=N-M) 
TZ VS TIME (GG@1000K) 
TIME( UN=SEC) , BX(UN=W-M) 
BX VS TIME (GG@1000K) 
TIME( UN=SEC) , BY(UN=W-M) 
BY VS TIME (GG@1000K) 
TIME( UN=SEC) , BZ(UN=W-M) 
BZ VS TIME (GG@1000K) 
TIME( UN=SEC) , PHID( UN=' RAD/SEC'’ ) 
PHID(RADSEC) VS TIME (GG@1000K) IX=95 IY=98 IZ=3. 00 
TIME( UN=SEC) , TD(UN=' RAD/SEC’ ) 
TD(RADSEC) VS TIME (GG@1000K) IX=95 IY=98 IZ=3.00 
TIME(UN=SEC) , PSID(UN='RAD/SEC' ) 
PSID(RADSEC) VS TIME (GG@1000K) IX=95 IY=98 I2Z=3. 00 
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